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Abstract. Based on sequences of three mitochondrial and two nuclear-encoded genes, we examine genetic variation in the 
leaf tailed geckos Uroplatus sameiti and U. sikorae, and morphological and chromatic characteristics of the genetic clusters 
identified. The mitochondrial phylogeny reveals a puzzling diversity of 16 deep lineages (4 in U. sameiti and 12 in U. siko- 
rae) differing by 2.9-9.9% uncorrected pairwise distance in a fragment of the 16S rRNA gene. Populations from Analalava 
in the North and Zahamena in the Northern Central East form two mitochondrial lineages clustering with U. sameiti but 
being deeply divergent (>8% 16S distance); however, no information on their morphology is available. In U. sikorae, the 
mitochondrial lineages identified form several major geographic clades, two of which (from the northernmost and south- 
ernmost populations, respectively) received substantial support in the phylogenetic analysis. No instance of sympatry of 
two or more mitochondrial lineages was observed, precluding an unambiguous assessment of species status under the bio- 
logical species criterion without experimental approaches or detailed hybrid zone analyses. In the fragment of the nuclear 
encoded gene SACS we observed haplotype sharing between species and mitochondrial lineages, while in the fragment of 
KIAA1239 no haplotype sharing was detected although neither species nor mitochondrial lineages formed coherent phylo- 
groups in the respective network. A screening of colour patterns from live photos, partly of the genotyped individuals, 
confirmed a large variation within species and populations, with a possible sexual dichromatism where a longitudinally 
striped phenotype is restricted to males. All individuals from populations of the U. sikorae clade from the Southern Central 
East and South East had an unpigmented oral mucosa just like U. sameiti, while all other U. sikorae populations are char- 
acterized by a black oral mucosa pigmentation. The extremely strong phylogeographic structure in the U. sikorae complex 
without obvious species-level divergences is unprecedented for large-sized squamates in Madagascar and calls for further 
taxonomic scrutiny using phylogenomic approaches; and it exemplifies how the loss of any major block of the remain- 
ing rainforests in Madagascar will inevitably lead to a substantial loss of genetic diversity - even if often intraspecific — in 
rainforest-specialized species. 
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Introduction though the first species of the genus was only scientifi- 


Leaf-tailed geckos of the genus Uroplatus DUMERIL, 1806 
are among the most prominent reptiles of Madagascar due 
to their bizarre appearance (GLAW & VENCES 2007, WOL- 
LENBERG et al. 2011). The first published reports of these 
strange creatures date back to 1658 (FLACOURT 1658), al- 


cally named 141 years later (SCHNEIDER, 1792). Curious- 
ly however, their taxonomy and geographic variation has 
only been intensively and comprehensively studied during 
the last twenty years, revealing that the genus is charac- 
terized by a much higher species diversity than previously 
assumed, with a high number of range-restricted species. 
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Insights into the complex phylogeography of Uroplatus sikorae and U. sameiti 


In the period between the first description in 1792 and 
1990, the species U. fimbriatus SCHNEIDER, 1792, U. linea- 
tus DUMERIL & BIBRON, 1836, U. ebenaui BOETTGER, 1879, 
U. alluaudi MOCQUARD, 1884, U. phantasticus BOULENGER, 
1888, U. guentheri MOCQUARD, 1908, and U. sikorae BOETT- 
GER, 1913 were scientifically named. Five species, U. henkeli, 
U. malahelo, U. malama, U. pietschmanni, and U. sameiti, 
were then discovered and scientifically named between 
1990-2004 (BÖHME & IBISCH 1990, NUSSBAUM & RAXWOR- 
THY 1994, 1995, BÖHLE & SCHÖNECKER 2003), followed by 
first molecular assessments that revealed widespread cryp- 
tic diversity and led to the definition of numerous candi- 
date species (GLAW et al. 2006, GREENBAUM et al. 2007, 
RAXworRTHY et al. 2008, RATSOAVINA et al. 2013). Subse- 
guently, eight new species of leaf-tailed geckos were for- 
mally described and named within this genus (U. finiavana 
RATSOAVINA et al. 20112, U. fiera RATSOAVINA et al. 2015, 
U. fotsivava RaTSOAVINA et al. 2017, U. kelirambo Ratso- 
AVINA et al. 2017, U. finaritra RATSOAVINA et al. 2019a, 
U. fetsy Ratsoavina et al. 2019b, U. fangorn RATSOAVINA 
et al. 2020a and U. fivehy Ratsoavina et al. 2020a). 

Taking into account the latest species descriptions 
(RATSOAVINA et al. 2020a), the genus Uroplatus current- 
ly contains 21 nominal species which can be grouped into 
five morphologically distinct species groups (RATSOAVINA 
et al. 2013). Over the last decade, taxonomic research focus 
has been on deciphering the species diversity within the 
small-sized, leaf-mimicking forms in the U. ebenaui group, 
although the second most species-rich group, the U. fim- 
briatus group, is also known to contain numerous genetic 
lineages of uncertain taxonomic status (RATSOAVINA et al. 
2013). Based on morphological similarity and phylogenet- 
ic relationships, the U. fimbriatus group can be subdivided 
into the U. fimbriatus complex (U. fimbriatus, U. gigante- 
us), the U. henkeli complex (U. henkeli and the candidate 
species U. henkeli [Ca11] - currently in the process of being 
described), and the U. sikorae complex (U. sameiti, U. siko- 
rae). Of these taxa, the genetic variation in the U. fimbria- 
tus complex has been analyzed by GEHRING et al. (2018), 
and the taxonomic revision of the U. henkeli complex is 
being completed as well (CG. Aw et al. submitted). Here, we 
focus on the U. sikorae complex where RATSOAVINA et al. 
(2013) documented, both within U. sameiti and U. sikorae, 
a large number of deep mitochondrial lineages whose mor- 
phology and relationships are largely unknown. 

Previous studies (GREENBAUM et al. 2007, RAXWORTHY 
et al. 2008, RATSOAVINA et al. 2013) based on multiple mo- 
lecular markers have provided clear evidence for the mono- 
phyly of the U. sikorae complex, with two major clades that 
correspond to U. sikorae and U. sameiti. The morphological 
diagnosis of these species, however, turned out to be con- 
tentious. Uroplatus sikorae was scientifically named and 
described by BOETTGER (1913), on the basis of its purport- 
ed smaller size compared to U. fimbriatus, as well as some 
scalation features which however were not considered to 
be diagnostic by subsequent authors. Prior to 1989, it was 
thus disputed whether U. sikorae represents a species sep- 
arate from U. fimbriatus. The systematic revision of Uro- 


platus by BAUER & RUSSELL (1989) resurrected U. sikorae 
and considered it a distinct species based on morphologi- 
cal characters, i.e., differences in dermal flaps and colora- 
tion, and sympatric occurrence with U. fimbriatus. Subse- 
quently, BÖHME & IBIscH (1990) found evidence for two 
separate taxa they considered as subspecies: U. sikorae siko- 
rae (type locality near Andrangoloaka) represented by the 
population from “Périnet” (=Andasibe, close to Mantadia- 
Analamazaotra National Park, geographical coordinates: 
-18.9333, 48.4167) and other localities mainly at mid-eleva- 
tions (roughly between 600 and 1,500 m a.s.l.), and U. siko- 
rae sameiti, with the type locality Nosy Boraha (= Sainte 
Marie, -16.9173, 49.87399) from mainly lowland localities. 
These two taxa were diagnosed mainly by the pigmentation 
of their oral mucosa: black in U. sikorae and unpigmented 
(pinkish in colour) in U. sameiti. RAXWORTHY et al. (2008) 
found substantial molecular divergence between these taxa 
and elevated them to species level, but RATSOAVINA et al. 
(2013) found evidence for a more complex situation, with 
numerous deeply divergent mitochondrial lineages both 
within sikorae and sameiti, and numerous populations in 
the genetic sikorae lineage lacking a black oral mucosa. 

In this study we aim to (i) more comprehensively assess 
the genetic diversity within the U. sikorae complex based 
on mitochondrial DNA, (ii) document the geographical 
distribution of the main mitochondrial lineages identified, 
as well as (iii) screen for possible chromatic and morpho- 
logical characters that would allow a distinction of the ge- 
netic groups. Our goal is to provide a detailed review of the 
U. sikorae complex as a baseline for future studies on their 
biogeography and systematics, and for improved conserva- 
tion management. 


Materials and methods 
Sampling for molecular analysis 


For molecular phylogenetics, we focused on amplifying 
and sequencing three mitochondrial gene fragments: 12S 
ribosomal RNA (12S rRNA or 12S), 16S ribosomal RNA 
(16S rRNA or 16S), and NADH dehydrogenase subunit 4 
(ND4) including partial tRNA stretches. We used sequenc- 
es compiled by Ratsoavina et al. (2013) and complement- 
ed these with new sequences available from GenBank, 
and further new sequences generated for this study. Since 
the 12S fragment was already available for many samples 
(RATSOAVINA et al. 2013) we attempted to complement the 
12S alignment for as many further samples as possible. As 
a specific challenge for this study, U. sikorae and U. sameiti 
contain a large number of deeply divergent mitochondrial 
lineages, and neither the short 12S fragment nor the ND4 
fragment alone were able to reliably reconstruct their rela- 
tionships (see RATSOAVINA et al. 2013); furthermore, many 
samples of U. sameiti/sikorae sequenced by other studies 
(e.g., RAXWORTHY et al. 2008, RATSOAVINA et al. 2013) 
were sequenced (and uploaded to GenBank) only for ei- 
ther 12S or ND4, and the respective tissue samples in many 
cases were not available to us. This unbalanced availabil- 
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ity of seguences led to numerous phylogenetic artefacts in 
exploratory analyses of the combined multigene data set, 
since representatives of some regional clades were present 
for only one inference of relationships. New seguencing 
was therefore directed at obtaining seguences of all three 
gene fragments for at least one sample per main mitochon- 
drial lineage. 


Molecular and phylogenetic methodology 


Total genomic DNA was extracted following a standard salt 
extraction protocol after proteinase K digestion (BRUFORD 
et al. 1992). Polymerase chain reaction with standard cy- 
cling protocols (BAUER et al. 2011) was carried out using the 
following primers: a fragment of the 12S ribosomal RNA 
gene (12S) was amplified with 12SAr-L 5 -AAACTGG- 
GATTAGATACCCCACTAT-3 and 12SBr-H 5 -GAGGGT- 
GACGGGCGGTGTGT-3 (Parumgi et al. 1991), a frag- 
ment of the 16S ribosomal RNA gene (16S) with 16SAr-L 
5-CGCCTGTTTATCAAAAACAT-3 and  16SBr-H 
5-CCGGTCTGAACTCAGATCACGT-3 (PALUMBI et al. 
1991), and a fragment of the mitochondrial gene for NADH 
dehydrogenase subunit 4 and following tRNAs (ND4) with 
ND4 5 -CACCTATGACTACCAAAAGCTCATGTAGAA- 
GC-3 and LeutRNA 5 -CATTACTTTTACTTGGATTTG- 
CACC-3 (ARÉVALO et al. 1994). In addition, we sequenced 
fragments of the nuclear-encoded genes for sacsin (SACS) 
and Leucine-rich repeat and WD repeat-containing pro- 
tein (KIAA1239) following primers and nested PCR ap- 
proach of SHEN et al. (2012). 

PCR products were purified with Exonuclease I and 
Shrimp Alkaline Phosphatase digestion, and the purified 
products along with sequencing primers were shipped 
to LGC Genomics (Berlin) for sequencing on automated 
capillary sequencing instruments, with a select number 
of samples processed and sequenced in-house at Omaha's 
Henry Doorly Zoo and Aquarium. Sequences were quality- 
checked and poor-quality stretches manually trimmed in 
CodonCode Aligner (Codon Code Corporation). We used 
MEGA7 (Kumar et al. 2016) for initial sequence align- 
ment, exploratory phylogenetic inference, and for calcu- 
lating uncorrected p-distances between sequences. Newly 
generated sequences were deposited in GenBank under 
the following accession numbers: 00302235-00302269, 
00303766-00303827, and O0318161-0 0318166. 

Seguences and associated metadata were curated in a 
Microsoft Excel spreadsheet, exported as tab-delimited 
text, and used as input for Concatenator (VENCES et al. 
20222) which is part of the iTaxoTools project (VENCES 
et al. 2021) where each gene fragment was aligned with 
MAFFT (KATOH & STANDLEY 2013) and concatenated into 
a Nexus-formatted file. For phylogenetic inference we per- 
formed a partitioned Bayesian analysis of the concatenat- 
ed mitochondrial fragments (2,225 bp) with MrBayes 3.2 
(RONQUIST et al. 2012), defining three partitions (12S + 16S, 
with a GTR4G model; ND4 third codon positions with a 
GTR+G model; and ND4 first and second codon positions 
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with a K2P+G model, based on model testing under the 
Bayesian Information Criterion in MEGA7), running 20 
million generations, sampling every 1,000 trees, and dis- 
carding the first 25% of sampled trees as burn-in. 

We inferred haplotypes of each of the nuclear DNA 
fragments with the PHASE algorithm (STEPHENS et al. 
2001) using the software DnaSP (Version 5.10.3; LIBRA- 
DO & ROZAS 2009). We reconstructed a Maximum Likeli- 
hood tree with Jukes-Cantor substitution model (chosen to 
avoid overparametrisation) in MEGA7 and then used this 
tree, separately for each fragment, to build a haplotype net- 
work by entering the tree together with the alignment in 
the software Haploviewer, written by G. B. Ewine (http:// 
www.cibiv.at/-greg/haploviewer), which implements the 
methodological approach of SALZBURGER et al. (2011). 


Morphological analyses 


For this study, we focused on specimens housed at the Zoo- 
logische Staatssammlung Miinchen (ZSM), collected over 
the last 20 years by ourselves and other colleagues on nu- 
merous expeditions. Field number abbreviations of these 
specimens and all others used in our phylogenetic trees re- 
fer to tissue or specimen numbers of P.-S. Gehring (PSG), 
R. A. Nussbaum (RAN), M. Pabijan (MPFC), A. Rakotoar- 
ison (AND), A. P. Raselimanana (APR), F. M. Ratsoavina 
(RATE, KAE KIAN, M, URANO, ZAH, ERT), C. J. Rax- 
worthy (RAX), M. D. Scherz (MSZC), M. Vences and/or F. 
Glaw (MVTIS, FG/MV, ZCMV, FGZC), and D. R. Vieites 
(DRV). 

For a selected number of specimens in the ZSM collec- 
tion we undertook a screening of possibly diagnostic mor- 
phological characters, including morphometric measure- 
ments, scale counts, and counts of the serrations in the lat- 
eral skin flaps. We focused on analyzing only those speci- 
mens in relatively good state of preservation, largely adults, 
which also had been included in the molecular analysis or 
from the same sites from which we had molecular data 
available. Unfortunately, this meant that only a limited 
number of voucher specimens of each genetic lineage was 
available for examination. We also qualitatively screened 
all individuals for additional characters such as skin tex- 
ture, tubercles, and additional scale counts, but did not 
further report on these given that no consistent differenc- 
es were observed. The selected measurements and counts 
were taken by student assistants, where each measurement 
was taken by the same person but in a haphazard order of 
individuals. For the measurements and counts taken, see 
Table 1. 


Analyses of colour variation 


To assess possible chromatic differences between indi- 
viduals or populations assigned to genealogical groups 
(identified above through phylogeographic analyses), es- 
pecially between mitochondrial lineages of U. sikorae and 
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U. sameiti, we compared colour patterns of the body, eye 
and buccal mucosa based on photographs of specimens 
from the entire distribution range of the two species, in or- 
der to follow an integrative taxonomic approach (PADIAL 
et al. 2010). All photographs had been taken of live geckos 
in the field during day and/or night. Our data did not al- 
low to systematically identify possible differences between 
day and night coloration, and we therefore do not further 
elaborate on these (except for iris coloration; see below); 
however, in general, no fundamental differences between 
day and night coloration were detected. Photographs were 
taken in localities that comprised intact or slightly to mod- 
erately disturbed rainforest and ranged in elevation from 
about 0-1,550 m a.s.l. (Supplementary Table S1). 


Results 
Molecular differentiation 


The molecular dataset of concatenated mitochondrial 
genes consists of 192 ingroup sequences for the three tar- 
geted mitochondrial DNA fragments (12S: 80 sequences; 
16S: 35 sequences; ND4: 77 sequences) for a total of 120 
individual ingroup samples, including 84 sequences new- 
ly determined for this study. The phylogenetic tree calcu- 
lated from the 2,225 bp of concatenated sequences (Fig. 1) 
divides the U. sikorae complex into two main monophyl- 
etic groups corresponding to U. sameiti and U. sikorae as 
currently understood, but found no significant support for 
these clades (Bayesian posterior probability PP = 0.60 for 
U. sikorae and < 0.5 for U. sameiti). The tree also reveals a 
puzzling diversity of deep mitochondrial lineages within 
these two taxa, especially within U. sikorae (lineages siki- 
sik12, Fig. 1) but also in U. sameiti (sam13-sam16; Fig. 1). 
All of these mitochondrial lineages are allopatrically dis- 
tributed (Fig. 2). Although various deeper nodes were 
poorly supported, some geographical groupings received 
high support: (i) Lineages siki-sik4, here together named 
as U. sikorae SCE/SE (PP = 0.96) encompassed all popu- 
lations from the Southern Central East and South East, 
ranging from Beampingaratra in the south-eastern Ano- 
sy Chain to Ranomafana National Park, and encompass- 
ing both low-elevation sites near sea level (e.g., Manombo) 
and mid-elevation sites (e.g., Ranomafana National Park). 
(ii) Lineages siks-sik8 were placed in an unsupported 
group here called U. sikorae NCE/NE (PP < 0.5) and in- 
cluded populations from mid-elevation localities, ranging 
from the Andasibe region in the Northern Central East to 
the Marojejy Massif in the North East, and reaching the 
limits of the Sambirano region at Ampotsidy. (iii) Lineages 
siko-siki2 formed a highly supported group here called 
U. sikorae N (PP = 1.0), encompassing populations from 
the northern massifs of Tsaratanana, Manongarivo, Sora- 
ta, and Montagne d'Ambre. (iv) Within U. sameiti, lineages 
sam13 and sam14 (grouped together with PP = 1.0), includ- 
ing specimens from the type locality Nosy Boraha, were 
found in primarily low-elevation sites along Madagascar’s 
east coast, with the two clades encompassing populations 
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from Marolambo to Betampona, and from Tampolo to 
Nosy Mangabe, respectively. (v) Finally, the enigmatic lin- 
eages samı5 and sam16 show a very deep molecular differ- 
entiation compared to typical U. sameiti. 

A comparison of elevational distribution of the lineages 
identified (Supplementary Table S1) revealed for locali- 
ties confirmed by molecular data an elevational range of 
U. sameiti from 4-774 m above sea level, and for U. sikorae 
from 279-1,550 m a.s.l., with the low elevation U. sikorae 
sites below 600 m all referring to the SCE/SE clade, outside 
of the latitudinal range of U. sameiti, all other localities of 
U. sikorae are at least at an elevation of 600 m, and mostly 
> 900 m a.s.l. 

Mitochondrial sequence divergences between lineages 
were high. In the 16S gene, uncorrected pairwise distances 
amounted to 6.3-11.4% between U. sameiti and U. sikorae; 
2.9-6.9% among lineages within U. sikorae SCE/SE, 4.5- 
8.4% among lineages within U. sikorae N, 3.3-7.9% among 
lineages within U. sikorae NCE/NE, and 4.3-9.6% among 
the three main geographical clades. The Zahamena lineage 
of U. sameiti differed by 8.6-9.9% from other U. sameiti. 

Sequences of the nuclear-encoded SACS (1068 bp) and 
KIAA1239 (872 bp) fragments were available each for 25 
individuals of the U. sikorae complex, covering two mi- 
tochondrial lineages of U. sameiti (sam13 and samı4) and 
seven lineages of U. sikorae (sik1, siks, sik7, sikg-siki2) that 
represent the main U. sikorae groups N, NCE/NE, and SCE/ 
SE. Although alleles (haplotypes) of the various mitochon- 
drial lineages showed some separation also in the nuclear 
gene haplotype networks (Fig. 3), several instances of al- 
lele sharing were apparent. Notably, in the SACS network, 
allele sharing was observed between U. sameiti from the 
more southern samı4 with an U. sikorae individual from 
sik7 (from Fierenana), and also other alleles from U. sikorae 
were placed in a separate phylogroup with U. sameiti indi- 
viduals, separated by a minimum of 5 mutational steps from 
other U. sikorae. In KIAA1239, no haplotype sharing be- 
tween any mitochondrial lineage was observed, but neither 
the two species nor any of the geographical clades formed 
coherent, separate phylogroups, and quite a large number 
of alleles were observed within lineages that differed by nu- 
merous mutational steps (up to 26 steps in sik12). 


Morphological characteristics 
of U. sameiti and U. sikorae 


Morphological data for representative individuals of 
U. sameiti and U. sikorae are given in Table 1. In the follow- 
ing account we compare these with general information 
summarized from the literature (BÖHME & IBISCH 1990, 
GLAW & VENCES 2007, GEHRING 2020). Snout-vent length 
(SVL) in U. sameiti is 103-110 mm in the three presum- 
ably adult females available to us, but is known to reach 
110-130 mm in this species. In U. sikorae, SVL ranges from 
97-120 mm in males and 93-116 mm in females in our sam- 
ples, but is known to reach up to 123 mm. Counting supra- 
labial and infralabial scales in these geckos is not easy be- 


Insights into the complex phylogeography of Uroplatus sikorae and U. sameiti 


URANO 4.13 
URANO 4.45 


URANO 4.26 
URANO 4.09 


0.02 MVTIS 31715 
PSG 2659 Ran 
1 IL FGMV 2002.31: 


KAF 137 Kianjavato 
0.99 KAF 187 Kianjavato 
KIAN 125 Kianjavato 
KIAN 144 Kianjavato 
1 KAF 125 Kianjavato ES 


Ee 


KL 
Le 
x~ 
S 
Q 
MPFC 508 Ran m 
Lë 
CH 
m 
~ 
Lë 
m 


KAF 047 Kianjavato 
1 KAF 338 Andalangina 
0.96 KAF 340 Andalangina 
KAF 265 Kianjavato 
1 |L_ KAF 153 Kianjavato 
KAF 154 Kianjavato 
AND 051 Manangotry 
AND 008 Manangotry 
AND 048 Manangotry 
MPFC 420 Beampingaratra 


0.98 


0.99 


= 


RAX 3644 Betaolona 
MSZC 117 Ampotsidy 
0.97 MSZC 131 Ampotsidy 
MSZC 121 Ampotsidy 
PSG 2523 Ambatovy 
1 PSG 2524 Ambatovy 
1 FGZC 4351 Anjozorobe 
FGZC 4386 Anjozorobe 
FGMV 2002.2232 Fierenana 
St FGMV 2002.951 Andasibe 
FGMV 2002.952 Andasibe 
APR 7636 Ambohitantely 


IN /3DN 2040yis N 


163 


PHILIP-SEBASTIAN GEHRING et al. 


0.64 


FGMV 2002.3011 
RAN 38129 Mont 


FRT 1009 M 


FGZC 4916 Cc 
a 
~ 
S 
Q 
Wl 
2 

0.63 
S 
v 
3 
PSG 183 Vohibola DD 
Gë 


PSG 203 Vohibola 
PSG 325 Sahafina 
PSG 1285 Ambodisavoka 
PSG 1379 Marolambo 
RAX 7700 Betampona 
PSG 324 Sahafina 
PSG 325 Sahafina 
PSG 1285 Ambodisavoka 
PSG 1379 Tanin Ny Betsileo 
MPFC .68 Camp Ambatofotsy 
MPFC 122 Vohimana sam14 


RAX 4483 Analalava 


Figure 1. Bayesian phylogenetic tree (5096 majority-rule consensus with all compatible bifurcations shown) based on DNA seguences 
of concatenated fragments of the mitochondrial genes for 12S rRNA, 16S rRNA, and ND4 of samples of the Uroplatus sikorae complex. 
Values at nodes indicate posterior probabilities (PP) > 0.95. 
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cause they are small, and become tiny towards the rictus. 
We here counted all scales around the upper and lower lip, 
including the very small ones in the corner of the mouth, 
and thus possibly obtained slightly higher counts than re- 
ported elsewhere in the literature. We observed 22-29 in- 
fralabials and 25-33 supralabials in U. sameiti, and 25-33 
infralabials and 29-39 supralabials in U. sikorae (Table 1). 
Overall, this might indicate a trend of more labial scales 
in U. sikorae, and within U. sikorae the highest values oc- 
curring in the northern mitochondrial lineages (Table 1). 
However, variation of this character was high across the 
specimens examined, and due to low sample sizes for most 
lineages we refrained from any statistical analysis. 


In individuals of both species, the body is dorsally flat- 
tened and the animals possess dermal fringes along their 
flanks, under their chin and on limbs, which clearly identi- 
fies them as members of the U. fimbriatus group. We report 
counts of the number of serrations of these fringes (Ta- 
ble 1) along head, flanks, and limbs, but did not encoun- 
ter any clear and consistent differences between species or 
mitochondrial lineages. The head is also strongly flattened 
and ends in an elongated snout. The rostral is not split. 
Scalation of the dorsal body surface is heterogenous and 
consists of more or less round granular scales that are in- 
terrupted by enlarged tubercle scales at irregular intervals, 
which are pointed or cone-like, and distinctly higher than 
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Figure 2. Geographical distribution of species and mitochondrial lineages in the Uroplatus sikorae complex, as assessed in the present 
work. Colours are according to the mitochondrial lineages in the phylogenetic tree (Fig. 1). 
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granular scales. The belly scales are granular, smooth, and 
relatively large, with a trend of being larger in the center 
of the belly. Transverse counts of belly scales at midbody 
ranged from 69-95, without a clear trend of differences be- 
tween species or mitochondrial lineages (Table 1). On the 
upper eyelid, there are 3-5 superciliary lobes. The third fin- 
ger bears 8-10 scansorial lamellae in both species (Table 1). 
At the tip of the tail there is no distinct notch and there are 
1-2 dermal spines at the root of the tail. The tail is elon- 
gated, dorsally flattened and two lateral dermal fringes run 
towards the rounded tip of the tail. 

Genital morphology was not assessed by us, but we 
here review available data from the literature. Hemipenis 


descriptions are available for both taxa in BOHME (1988) 
and BOHME & IBISCH (1990) from Andasibe and Nosy Bo- 
raha, respectively. From these descriptions, hemipenes 
both of U. sameiti and U. sikorae appear to be character- 
ized by a shield-like bulge of the upper border of the asul- 
cal side of the hemipenis that exceeds the serrated apical 
lobes. This bulge is covered by deep, honeycomb-like caly- 
ces on the asulcal side and embraces the truncus left and 
right in the upper third of the sulcal side. The presence of 
this bulge differentiates the male genitals of the U. sikorae 
complex from those of U. fimbriatus and U. henkeli, but no 
differences between U. sameiti and U. sikorae were noted 
(BOHME & IBISCH 1990). 
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Figure 3. Haplotype networks reconstructed from phased alleles of fragments of the nuclear-encoded SACS (1068 bp) and KIAA1239 
(872 bp) genes for 25 individuals of the Uroplatus sikorae complex (thus, each individual is represented with two sequences in the 
analysis of each gene). Haplotype colours and names correspond to mitochondrial lineages in Figure 1. 
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Comparison of colour patterns 


Photographs of 73 living individuals of U. sikorae (44 fe- 
males, 19 males, 10 unknown sex) from Andrafainkoa, Am- 
bolopatrika, Ampotsidy, Anjozorobe, Andasibe (Anala- 
mazaotra), Beampingarata, Kianjavato, Makira, Marojejy, 
Marolambo, Maromizaha, Montagne d'Ambre, Ranomafa- 
na, Sorata and Tsararano were available for analysis. Thus, 
all mitochondrial lineages of U. sikorae are photographical- 
ly covered by representative individuals, and many of these 
corresponded to the specimens identified via molecular 
data. For U. sameiti photographs of 26 individuals (14 fe- 
males, 14 males and 2 of unknown sex) from Ambodiriana, 
Betampona, Masoala, Nosy Boraha, Nosy Mangabe, Sa- 
hafina, Tampolo, Vohibola and Vohimana were analyzed. 
All photographed individuals belonged to lineages sam13 
and sam14, unfortunately no pictures of animals from Za- 
hamena (lineage sam1s) or Analalava (sam16) were avail- 
able to us. Many of the photographed individuals corre- 
sponded to the specimens identified via molecular data. 
Additional photographs of U. sameiti and U. sikorae, which 
in general comply with the pattern described here, can be 
found in the publications of SVATEK & VAN DUIN (2002), 
SCHONECKER (2008) and GEHRING (2020) as well as on the 
internet (e.g. iNaturalist). 

Specimens of U. sikorae exhibit the following dorsal col- 
oration (an exemplified overview is given in Figs 4-8): the 
dorsal coloration is highly variable across individuals and 
mitochondrial lineages; it generally mimics moss- and li- 
chen-covered trunks of rainforest trees. There are mossy- 
green, lichen- or bark-coloured, reddish or almost com- 
pletely white specimens. There can be blotches, spots or 
stripes present as patterns. Longitudinal stripes on the 
back, neck and head have only been observed in males 
which may represent a difference between sexes; however, 
males can also be without a striped pattern (e.g. Figs 5C, 
5E). We were unable to detect any consistent differences 
in dorsal coloration between the mitochondrial lineages of 
U. sikorae that would allow unambiguous assignment. In 
animals from the northern populations (Montagne dAm- 
bre, Tsaratanana; lineages siko and siki1) a fine reticulation 
is common that we have not so far observed in individuals 
of the populations from the southeast coast or the central 
highlands (e.g. Figs 8A, 8E, 8I, 8J). 

The ventral surface of all populations is light grey to 
white, and single, irregularly dispersed dark spots can be 
found from the tip of the snout to the tip of the tail. Ven- 
tral dermal fringes of the tail are densely covered with dark 
spots, so they form larger blotches. 

As already shown in Ratsoavina et al. (2013), the col- 
oration of the oral mucosa is not constant within U. sikorae 
and differs between main geographical clades. In all south- 
ern populations (lineages siki-sik4: Fig. 1a), the oral mu- 
cosa is unpigmented (Figs 4 and 5). In the clades from the 
central highlands and northern Madagascar (lineages siks— 
sik12, Fig. 1a and 1b), the oral mucosa is darkly pigmented 
(Figs 6-8). The tongue is strikingly yellow coloured in two 
photographed individuals from Ampotsidy (sik6; Fig. 6C). 


According to the photographs examined, U. sameiti 
exhibits the following dorsal coloration (Fig. 9): the dor- 
sal coloration is highly variable and generally mimics tree 
trunks of rainforests covered in mosses and lichens. There 
are mossy-green, lichen- or bark-coloured, greyish or al- 
most completely white-coloured specimens. Individuals 
show a variety of different patterns such as a fine reticu- 
lation, dots or even stripes. Also in this species, a longi- 
tudinally striped pattern was only observed in males. The 
ventral surface of all populations is light grey to white and 
single, irregularly dispersed dark spots can be found from 
the tip of the snout to the tip of the tail. Ventral dermal 
fringes of the tail are densely covered with dark spots, so 
they form larger blotches. No clade-diagnostic colour pat- 
terns were observed between lineages sami3 and sami4. 
The oral mucosa is constantly unpigmented in all popula- 
tions of U. sameiti. 

In both species, the eyeballs are encircled by a bright 
yellow dermal ring that becomes visible when the animals 
perform their defensive behaviour and protrude their eyes 
from the sockets. The outer iris area is yellow to light grey 
or reddish, sometimes with a bluish shade, by day, while at 
night it is yellow at the outer border and reticulated with an 
irregular network of fine veins that are condensed towards 
the pupil, giving the iris a brown colour. 

To summarize, an assignment to genetically defined 
species or mitochondrial lineages in the U. sikorae com- 
plex is not possible based on coloration. The distinction 
between U. sameiti and U. sikorae based on the pigmenta- 
tion of the oral mucosa is not unambiguous since an un- 
pigmented oral mucosa is also found in the SCE/SE mi- 
tochondrial lineages of U. sikorae. Northern and southern 
U. sikorae lineages can thus be distinguished based on the 
coloration of the oral mucosa (black in the north, pink in 
the south). Note, however, that support for these clades is 
only tentative, and this feature may not therefore be val- 
uable as a trait for species delimitation. An ungenotyped 
specimen with a pink oral mucosa cannot be reliably as- 
signed to either the southern U. sikorae or U. sameiti. 


Discussion 


In this study we corroborate previous, preliminary data 
(e.g., RATSOAVINA et al. 2013) of strong phylogeographic 
structure in the Uroplatus sikorae complex. Because our 
molecular sampling of the three mitochondrial gene frag- 
ments analyzed (12S, 16S, ND4) were unevenly distribut- 
ed despite additional sequencing efforts, our data set was 
unsuitable for an objective lineage-delimiting tool such as 
ASAP (PUILLANDRE et al. 2021); we are, however, confident 
that the majority of ad-hoc mitochondrial lineages identi- 
fied (with 16S distances of 2.9-9.9% among them) would 
also be identified by objective criteria as putative species- 
level units. 

Recent phylogeographic studies have uncovered high 
genetic variation in numerous widespread reptile spe- 
cies of Madagascar (BOUMANS et al. 2007, FLORIO et al. 
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Figure 4. Representative specimens and colour variation of U. sikorae of the Ranomafana population (geographical clade U. sikorae 
SCE/SE; mitochondrial lineage sik1). (F) corresponds to sample PSG 2654. All photos by the authors. 
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Figure 5. Representative specimens and colour variation of Uroplatus sikorae (clade U. sikorae SCE/SE) of the populations of (A-D) Ki- 
anjavato (mitochondrial lineage sik2), (E-G) Beampingarata (lineage sik3), and (H-I) Andohahela populations. (F) corresponds to 
sample MPFC 420 in the phylogenetic tree. All photos by the authors. 
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Figure 6. Representative specimens and colour variation of Uroplatus sikorae of the mitochondrial clade U. sikorae NCE/NE. (A, B) In- 
dividuals of the Makira population (mitochondrial lineage sik5); (C-E) individuals from Ampotsidy (sik6); (F-H) individuals of the 
Marojejy population (sik5). (I, J) individual from Tsararano; (K, L) individuals of the Ambolopatrika population. Tsararano and Am- 
bolopatrika are located between the Anjanaharibe-Sud and Marojejy massifs. Photos I-L: F. ANDREONE; other photos by the authors. 
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Figure 7. Representative specimens and colour variation of Uroplatus sikorae of the mitochondrial clade U. sikorae NCE/NE, lineage 
sik7. (A, B, D, E, G) Individuals of the Andasibe-Analamazaotra population; (C) female from Anjozorobe; (F) female of the Maro- 
mizaha population, a rainforest fragment close to Andasibe. Photos A, G: A. HARTIG; F: F. ANDREONE; other photos by the authors. 
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Figure 8. Representative specimens and colour variation of Uroplatus sikorae of the mitochondrial clade U. sikorae N. (A-G) Indi- 
viduals of the Montagne d'Ambre population (sik9); (H-J) individuals of the Tsaratanana population (sik11); (K, L) individuals from 
Sorata; (M) individual from Andrafainkona. Photos B, E: A. Hanne: C: T. ALTHAUS; tb ANDREONE; other photos by the authors. 
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Figure 9. Representative specimens and colour variation of Uroplatus sameiti. (A, B, D) Individuals from Vohibola; (C, E) individuals 
of the Sahafina population (mitochondrial lineage sam14); (F, J) individual from Nosy Mangabe; (G-I) individuals of the population 
from the type locality Nosy Boraha (sam13); (K) individual from Ambodiriana (sam13). All photos by the authors. 
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2012, VENCES et al. 2014, GRBIC et al. 2015). In Uroplatus, 
the majority of recent molecular studies dealt with small- 
sized species which often are range-restricted; however, 
the large-sized U. giganteus and, to a lesser degree, U. fim- 
briatus do show substantial mitochondrial differentiation 
across their range (GEHRING et al. 2018). The two wide- 
spread species of small-sized leaf-tailed geckos (U. ebenaui 
and U. phantasticus) are equally subdivided into various 
deeply divergent lineages, at least from a mitochondri- 
al perspective (RATSOAVINA et al. 2012, 2019b). Still, the 
number of 12 deep mitochondrial lineages in U. sikorae 
represents an exceptional amount of genetic variation 
even in Madagascar’s herpetofauna, and it is not straight- 
forward to find an explanation for this pattern. Among 
possible factors, it would be worthwhile to investigate in 
depth (i) the detailed pattern of exceptionally high mito- 
chondrial substitution rates in geckos, and in particular 
in Uroplatus, i.e., by quantifying these substitution rates 
at different genes and codon positions as well as satura- 
tion levels of mitochondrial divergences in comparison to 
other animals; (ii) the possibility of low effective popula- 
tion sizes in Uroplatus, which in combination with repeat- 
ed bottlenecks could lead to high substitution rates and 
high phylogeographic fragmentation; (iii) the possibility 
of small individual home ranges and high site fidelity of 
Uroplatus which also could contribute to phylogeograph- 
ic fragmentation; and (iv) possible range fragmentations 
and habitat contractions of Uroplatus during past episodes 
of climate change, given that these geckos are strictly de- 
pendent on rainforest which may have experienced phas- 
es of contraction over the Pleistocene (e.g., WILME et al. 
2006). 

Range sizes in Madagascar’s reptiles decrease with body 
size (BROWN et al. 2016), and the microendemism observed 
in the Uroplatus ebenaui group has been hypothesized to 
be related to their smaller size (RATSOAVINA et al. 20204). 
Such an explanation cannot be put forward for the U. siko- 
rae complex, which is composed of relatively large-sized 
geckos, and alternative explanations are therefore required 
to understand their strong phylogeographic structure. One 
key might be their elevational occurrence. In fact, from a 
visual exploration of the distributional information, it ap- 
pears that the lowland lineages sam13 and sam14, belong- 
ing to U. sameiti, have colonized a wider area of Madagas- 
car’s eastern coast than any of the U. sikorae lineages that 
all occur at higher elevation. Also the two species of the 
U. fimbriatus complex, with a less pronounced phylogeo- 
graphic structure (RAXWORTHY et al. 2008, GEHRING et 
al. 2018), occur mostly at lower elevation. If such a trend 
would be confirmed by a meta-analysis relating range size 
to elevation in Madagascar’s rainforest biota, it could be an 
indication that low-elevation rainforests would have been 
less impacted by long-term fragmentation during past cli- 
mate fluctuations than those at mid- and high-elevations 
- despite such a scenario being counter-intuitive in light of 
the riverine barrier mechanism which postulates that riv- 
ers constitute stronger barriers to gene flow in the lowlands 
where they are wider (e.g., GEHRING et al. 2012). 
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The mitochondrial phylogenetic analysis itself is ham- 
pered by the large inventory of distinct lineages in the 
U. sikorae complex (which might not even be complete yet, 
as several areas have not been intensively sampled). There- 
fore, it is challenging to subdivide the puzzling diversity 
of lineages into larger clades and infer their evolutionary 
and biogeographic history. Our efforts to complement the 
12S and ND4 sampling of Ratsoavina et al. (2013) with a 
third gene fragment (16S) for representatives of almost all 
lineages led to an improvement of phylogenetic resolution, 
suggesting that an extended gene sampling might be able 
to fully resolve the phylogenetic tree of these geckos — at 
least from a mitochondrial perspective. Our tree (Fig. 1) 
resolves two main clades of U. sikorae with substantial sup- 
port, one that groups the southernmost lineages (the SCE/ 
SE clade; with unpigmented oral mucosa) and one that 
groups the northernmost lineages (the N clade), support- 
ing a pattern of geographical differentiation. However, the 
third U. sikorae clade (NCE/NE) is not adequately support- 
ed in our analysis. It is obvious that the available data are 
insufficient to conclusively fully resolve the relationships 
among lineages in the U. sikorae complex. Furthermore it 
is uncertain to which extent the mitochondrial signal may 
be blurred by introgressive hybridization or the presence 
of “ghost lineages” differentiated in their mitochondrial 
but not nuclear genome. Phylogenomic and population 
genomic data set will be necessary both for improved phy- 
logenetic resolution and species delimitation, especially 
by targeted sampling of contact/hybrid zones among line- 
ages to understand the extent of genome-wide admixture 
among them (e.g., DUFRESNES et al. 2021). 

Do any of the lineages or major geographical clades 
of U. sikorae or U. sameiti represent distinct, scientifical- 
ly unnamed species? Our data do not provide unambig- 
uous support for this hypothesis. In no case did we find 
evidence for two mitochondrial lineages occurring syn- 
topically which — in combination with lack of admixture 
in nuclear-encoded genes or morphological differences - 
could provide strong evidence for the existence of at least 
two separated evolutionary lineages using the biological 
species criterion. The two nuclear genes analyzed, further- 
more, did not provide a fully concordant signal with the 
mitochondrial tree, and neither the two currently accept- 
ed species (U. sameiti and U. sikorae) nor any of the ma- 
jor geographical clades belonged to separate phylogroups 
without haplotype sharing in the SACS network (Fig. 2). 
In the KIAA1239 network, no haplotype sharing among 
lineages was observed, but again the haplotypes did not 
form coherent phylogroups and the number of individu- 
als sequenced per lineage was relatively low. In addition, 
we did not find any morphological character unambigu- 
ously diagnostic for any species, geographical clade, or lin- 
eage; the sole informative character identified, black pig- 
mentation of the oral mucosa, characterizes specimens of 
U. sikorae from the NCE/NE and N clades which may not 
be each other’s closest relatives according to the mitochon- 
drial tree; and the absence of pigmentation is shared by 
U. sameiti and U. sikorae SCE/SE. 
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Our data are of importance to guide future studies of 
morphological differentiation in the U. sikorae complex 
which should focus on identifying combinations of mor- 
phological and perhaps osteological characters to differ- 
entiate the main genetic lineages. Importantly, individuals 
from different geographical areas should not be uncriti- 
cally combined in any comparison, and a reference to the 
mitochondrial clades identified herein should be included. 
From our screening and the data presented by BOHME & 
IBISCH (1990), it is unlikely that diagnostic characters in ex- 
ternal colour pattern exist, and the probable sexual dimor- 
phism in the occurrence of longitudinal stripes is probably 
a common pattern of various large-sized leaf-tailed geckos. 
This has also been documented in U. henkeli where this pat- 
tern occurred in males only (FOLEY & PFAFF 2005). Con- 
sidering our experience with other morphologically cryp- 
tic geckos in Madagascar (e.g., Lygodactylus; Vences et al. 
2022b), counts of the tubercular scales along the body axis 
(dorsal and ventral) may be informative to distinguish be- 
tween species, but their counting is very time-consuming. 
Large-sized Uroplatus are the squamates — and possibly the 
amniotes — with the largest number of teeth (BAUER & Rus- 
SELL 1989), and we anticipate that tooth number might be 
a taxonomically valuable character also within the U. siko- 
rae complex. Of particular interest is the hemipenial struc- 
ture which - providing information on sexual selection - is 
a prime taxonomic character in squamates (BOHME 1988, 
PapDiaL et al. 2010). A more comprehensive analysis of sex- 
ually mature male specimens of the U. sikorae complex, of 
reliably known geographical provenance genetically as- 
signed to lineages, and with fully everted hemipenes (cur- 
rently not available to us), is needed to assess whether dif- 
ferences in genital structures — even if subtle - may exist 
among lineages and species of this complex. 

Two mitochondrial lineages within U. sameiti are par- 
ticularly enigmatic and worthy of future taxonomic revi- 
sion: sam15 and sam16. Only a single sample of sam16 from 
Analalava was sequenced by Raxwortuy et al. (2008) for 
which neither 16S or ND4 sequences nor any morphologi- 
cal information is available. For sam15, four samples from 
Zahamena were sequenced for ND4 by Ratsoavina et al. 
(2013) and, independently, one sample from the same lo- 
cality for 12S and other genes by RAxworruy et al. (2008). 
We here added ND4 for two additional samples, and 12S 
and 16S sequences for several of our Zahamena samples, 
which allowed us to confirm that the sequences of Rax- 
WORTHY et al. (2008) and RaTSOAVINA et al. (2013) refer to 
the same mitochondrial lineage. This population is geneti- 
cally highly divergent (16S distance >8% to all other sam- 
ples in the U. sikorae complex) and its placement sister to 
U. sameiti sam13/sam14 received no statistical support. This 
deep lineage (possibly along with the one from Analalava 
in the north) may well represent a new species of leaf-tailed 
gecko, but since we had neither life photos nor voucher 
specimens available for examination, its status cannot be 
reliably assessed. Strikingly, also in the U. ebenaui group of 
small-sized leaf tail geckos, Zahamena harbors an appar- 
ently microendemic and genetically highly distinct candi- 


date species (U. ebenaui (Carol). In general, the herpeto- 
fauna of Zahamena National Park (GooDMAN et al. 2018) 
is among the taxonomically least explored amphibian and 
reptile communities of Madagascar and definitely requires 
future inventory work, of leaf tail geckos as well as other 
taxa. 

It is striking how poorly known the taxonomy as well as 
the natural history of the large-sized species of leaf tailed 
geckos is, considering the high value of these animals as 
flagship species for ecotourism and nature conservation in 
Madagascar (e.g., WOLLENBERG et al. 2011): Most data on 
the natural history and reproduction of Uroplatus available 
to date derive from captive observations (as summarized 
by GEHRING 2020), and to our knowledge no dedicated 
long term field study focusing on these aspects has been 
carried out on any Uroplatus. Despite the past and present 
importance of Uroplatus in the pet trade (e.g., RAXWORTHY 
& NUussBAUM 2000, ToDD 2011), not even thorough popu- 
lation density data based on mark-recapture experiments 
have been published for most species (except transect den- 
sity estimates for U. giganteus: INGADY 2011). 

Due to this complex taxonomic situation within the 
U. sikorae complex, information regarding the conserva- 
tion status of U. sameiti and U. sikorae must be considered 
preliminary and should be urgently re-assessed within the 
IUCN Red List framework. The current distribution range 
for U. sameiti for instance is given with an area of 52,955 km? 
of lowland areas along Madagascars east coast, therefore, 
the species is currently listed as “Least Concern” by the 
IUCN (RatTsoavina et al. 2020b). However, the assess- 
ment also lists all distribution points of the U. sikorae SCE/S 
clades, since genetic information was not yet available at the 
time of the assessment and all populations with an unpig- 
mented oral mucosa were identified as U. sameiti. Based on 
our new findings, it must be assumed that the actual size 
of the assumed distribution range for U. sameiti is there- 
fore only about half as large, and especially the lowland and 
coastal forests of Madagascar’s east coast that are under im- 
mense pressure of deforestation and fragmentation of habi- 
tat. Nevertheless, U. sameiti still occurs in numerous pro- 
tected areas (e.g. Vohimana, Betampona, Nosy Mangabe, 
Ambodiriana, Masoala), so it can be assumed that there is 
no immediate threat of extinction to the species. 

Currently, a distribution area of 77,103 km? is given for 
U. sikorae in the Red List database and the species is there- 
fore considered as “Least Concern” by the IUCN (Ratso- 
AVINA et al. 2011b). Recorded localities of the species as 
currently defined include some protected areas such as 
Analamazaotra, Mantadia, Maromizaha and Anjozorobe. 
Again, the problem is that the assessment arose before 
the remarkable genetic diversity in the complex was rec- 
ognized, and it is currently still questionable whether the 
genetic differentiation will lead to further splitting of the 
species after further investigation. If, after thorough taxo- 
nomic reassessment, U. sikorae sensu lato is split into mul- 
tiple species, these will necessarily have smaller ranges 
than the complex as a whole, and may therefore be threat- 
ened under the range-related Criterion B of the IUCN Red 
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List Categories and Criteria (IUCN, 2012). Treatment of a 
species complex with multiple, known potential species- 
level lineages as Least Concern could result in extinction 
of species going unnoticed, and on this grounds it has been 
argued that such species complexes should be assessed as 
Data Deficient, with emphasis on the need for a compre- 
hensive taxonomic revision before an accurate threat sta- 
tus assessment can be made (SCHERzZ et al. 2019). However, 
due to the puzzling taxonomic situation of U. sikorae de- 
spite the extensive data presented in this paper, we here re- 
frain from an IUCN threat status suggestion. 

The correct identification of species to be exported by 
the authorities’ officials for the pet trade is another prob- 
lem. Between 2005 and 2019, annually 1,500-2,000 in- 
dividuals of U. sikorae were released for export (UNEP- 
WCM 2015) and 1,391 individuals were exported between 
2019 and 2022 (October) (CITES Trade Database 2022). 
In the years 2014 and 2015, 500 individuals of U. sameiti 
were respectively released for export by Malagasy author- 
ities (UNEP-WCM 2015). Between 2019 and 2022 (Octo- 
ber) 498 individuals were exported (CITES Trade Database 
2022). In many cases it is hardly possible for the authorities 
and breeders to discriminate between U. sikorae and its sis- 
ter taxon U. sameiti, and differentiating among genetic lin- 
eages is completely impossible without precise information 
on collection site; avoiding unintentional cross-breeding 
among these lineages is therefore a big challenge. 

According to the data presented herein, populations of 
the U. sikorae complex from almost every remaining rain- 
forest block in Madagascar show substantial genetic diver- 
gence. Even if this may not be indicative of cryptic taxo- 
nomically-relevant diversity, it still exemplifies how further 
loss of primary rainforest habitat in Madagascar will inevi- 
tably lead to a loss of hitherto unknown and undescribed 
species and intraspecific genetic diversity, and calls for cau- 
tion when assembling captive populations for conservation 
breeding from specimens of unknown provenance. 


Acknowledgments 


We are grateful to numerous colleagues, students and guides for 
companionship and help during numerous field expeditions in 
Madagascar, carried out between 2001-2018, in particular GEN- 
NARO APREA, LAWRENCE BALL, PARFAIT BORA, JAMES BORRELL, 
HILDEGARD ENTING, OLIVER HAWLITSCHEK, ANGELIKA KNOLL, 
JORN KOHLER, ELLA Z. LATTENKAMP, NICOLA LUTZMANN, KON- 
RAD MEBERT, DENISE H NOMENJANAHARY, MACIEJ PABIJAN, 
DUNCAN PARKER, JIM and CAROL PATTON, DAVID PROTZEL, 
MARTA PUENTE, JEANNENEY RABEARIVONY, LILIANE RAHARI- 
VOLOLONIAINA, E. RAJERIARISON, T. RAJOAFIARISON, MARIUS 
RAKOTONDRATSIMA, ANDOLALAO RAKOTOARISON, SAFIDY M. 
RAKOTOMALALA, ONJA RANDRIAMALALA, LALAINA RANDRIA- 
MANANA, F. RANDRIANASOLO, ROGER-DANIEL RANDRIANIAINA, 
S. RASAMISON, ANGELUC RAZAFIMANANTSOA, ELIDIOT RAZAFI- 
MANDIMBY, JARY H. RAZAFINDRAIBE, THOMAS STARNES, RICKY 
TIAVINA, MEIKE TESCHKE and Davip R. ViEITEs. We also wish 
to thank THOMAS ALTHAUS, YVONNE ALTHAUS, FRANCO AN- 
DREONE, HANS-PETER BERGHOF, ARNE HARTIG, DAVID PROT- 
ZEL, GONÇALO M. Rosa and SEBASTIAN WOLFE who provid- 


176 


ed comparative photographs. Thanks are also due to the BDo8 
“Vertebrate Morphology” course of 2018 at TU Braunschweig 
(Kim-LarissA BAARS, JACQUELINE BARNBECK, JOHANNES FITZ- 
KE, JESSICA GEISSLER, LORENZ HELDT, SABRINA KABOTH, LEN- 
NART KRAUSE, CHRISTIAN LASCHEIT, MARCO LEITERHOLT, ANNE 
SCHELER, DANIELA SCHOOP, MELANIE SCHULZ, Nico WESCHE) 
who painstakingly contributed to performing scale counts and 
measurements under the supervision of MV and MDS. Thanks 
to JOSEF SCHMIDTLER for providing information on the works 
of J.G. SCHNEIDER. Research was carried out in the framework 
of collaboration accords between the Zoological Institut of TU 
Braunschweig, the Zoologische Staatssammlung Miinchen, the 
Mention Zoologie et Biodiversité Animale of the Université 
dAntananarivo, and the Ministere de PEnvironnement, des Eaux 
et des Foréts of the Republic of Madagascar. We are grateful to 
the Malagasy authorities and to Madagascar National Parks for 
research, collection, and export permits. 


References 


AREVALO, E., S. K. Davis & J. W. Sites JR. (1994): Mitochondri- 
al DNA-sequence divergence and phylogenetic relationships 
among eight chromosome races of the Sceloporus grammicus 
complex (Phrynosomatidae) in Central Mexico. - Systematic 
Biology, 43: 387-418. 

Bauer, A. M. & A. P. RUSSELL (1989): A systematic review of the 
genus Uroplatus (Reptilia: Gekkonidae) with comments on its 
biology. - Journal of Natural History, 23: 169-203. 


BAUER, A., F. Giaw, P.-S. GEHRING & M. Vences (2011): New 
species of Blaesodactylus (Squamata: Gekkonidae) from An- 
karafantsika National Park in north-western Madagascar. — 
Zootaxa, 2942: 57-68. 


BOHLE, A. & P. SCHONECKER (2003): Eine neue Art der Gattung 
Uroplatus Dumtril, 1805 aus Ost-Madagaskar (Reptilia: Squa- 
mata: Gekkonidae). - Salamandra, 39: 129-138. 


BÖHME, W. (1988): Zur Genitalmorphologie der Sauria: funktio- 
nelle und stammesgeschichtliche Aspekte. - Bonner Zoologi- 
sche Monografien, 27: 1-176. 


BOuME, W. & P. IBIscH (1990): Studien an Uroplatus. I. Der Uro- 
platus fimbriatus-Komplex. - Salamandra, 26(4): 246-259. 


BOETTGER, O. (1913): Reptilien und Amphibien von Madagascar, 
den Inseln und dem Festlande Ostafrikas (Sammlung Voeltz- 
kow 1889-1895 und 1903-1905). — In: VOELTzKOw A. (ed.): 
Reise in Ostafrika in den Jahren 1903-1905. Wissenschaftliche 
Ergebnisse. Dritter Band. Systematische Arbeiten. - Schwei- 
zerbartsche Verlagsbuchhandlung, Stuttgart. 


BOUMANS, L., D. R. VIEITES, F GLAw & M. VENCEs (2007): Geo- 
graphical patterns of deep mitochondrial differentiation in 
widespread Malagasy reptiles. - Molecular Phylogenetics and 
Evolution, 45: 822-839. 


Brown J. L., N. SILLERO, F. GLaw, P. Bora, D. R. VIEITES & M. 
Vences (2016): Spatial biodiversity patterns of Madagascar’s 
amphibians and reptiles. - PloS ONE, 11: e0144076. 

BRUFORD, M. W., O. HANOTTE, J. F. Y. BROOKFIELD & T. BURKE 
(1992): Single-locus and multilocus DNA fingerprint. - pp: 
225-270 in: HOELZEL, A. R. (ed.): Molecular Genetic Analy- 
sis of Populations: A Practical Approach. - IRL Press, Oxford. 


CITES Trade Database (2022): Compiled by UNEP-WCMC 


for the CITES Secretariat. - URL: trade.cites.org, accessed 
09.10.2022. 


Insights into the complex phylogeography of Uroplatus sikorae and U. sameiti 


DUEFRESNES, C., A. BRELSFORD, D. L. JEFFRIES, G. MAZEPA, T. SU- 
CHAN, D. CANESTRELLI, A. NICIEZA, L. FUMAGALLI, S. DUBEY, 
I. MARTINEZ-SOLANO, S. N. LITVINCHUK, M. VENCES, N. PER- 
RIN & P.-A. CROCHET (2021): Mass of genes rather than master 
genes underlie the genomic architecture of amphibian specia- 
tion. - Proceedings of the National Academy of Sciences of the 
U.S.A., 118: €2103963118. 


FLAcouRrg, E. DE (1658): Histoire de la Grande Isle de Madagas- 
car, composée par le Sieur de Flacourt, Directeur Général de 
la Compagnie Francaise de Orient et Commandant pour su 
Majesté dans la Dite Isle et Isles Adjacentes. - Alexandre Les- 
selin, Paris. 


FLORIO, A. M., C. M. INGRAM, H. A. RAKOTONDRAVONY, E. E. 
Louis & C. J. RaxworTHY (2012): Detecting cryptic specia- 
tion in the widespread and morphologically conservative car- 
pet chameleon (Furcifer lateralis) of Madagascar. - Journal of 
Evolutionary Biology, 25: 1399-1414. 

FOLey, S. C. & S. PFAFF (2005): Evidence of sexual dimorphism in 
neonate Henkel’s Leaf-tailed Gecko, Uroplatus henkeli. - Her- 
petological Review, 36: 156-157. 

GEHRING, P.-S. (2020): Leaf-tailed Geckos - The Complete Uro- 
platus. - Chimaira, Frankfurt am Main. 


GEHRING, P. S., M. PABIJAN, J. E. RANDRIANIRINA, F. GLAW & M. 
Vences (2012): The influence of riverine barriers on phylogeo- 
graphic patterns of Malagasy reed frogs (Heterixalus). - Mo- 
lecular Phylogenetics and Evolution, 64: 618-632. 


GEHRING, P. S., S. SIARABI, M. D. SCHERZ, F. M. RATSOAVINA, A. 
RAKOTOARISON, F. GLAW & M. Vences (2018): Genetic dif- 
ferentiation and species status of the large-bodied leaf-tailed 
geckos Uroplatus fimbriatus and U. giganteus. - Salamandra, 
54: 132-146. 

GLaw, F., J. Kosucu, F-W. HENKEL, P. SouND & W. BÖHME 
(2006): Genetic and morphological variation of the leaf-tailed 
gecko Uroplatus fimbriatus from Madagascar, with description 
of a new giant species. - Salamandra, 42: 129-144. 


Graw, F. & M. Vences (2007): A field guide to the amphibians 
and reptiles of Madagascar. Third edition. - Vences and Glaw 
Verlag, Cologne. 


GooDMAN, S. M., M. J. RAHERILALAO & S. WOHLHAUSER (eds) 
(2018): Les aires protégées terrestres de Madagascar: leur his- 
toire, description et biota. - Association Vahatra, Antanana- 
rivo. 


GRBIC, D., S. V. SAENKO, T. M. RANDRIAMORIA, A. DEBRY, A. P. 
RASELIMANANA & M. C. MILINKOVITCH (2015): Phylogeogra- 
phy and support vector machine classification of colour vari- 
ation in panther chameleons. - Molecular Ecology, 24: 3455- 
3466. 


GREENBAUM, E., A. M. BAugr, T. R. JACKMAN, M. VENCES & 
F. GLaw (2007): A phylogeny of the enigmatic Madagascan 
geckos of the genus Uroplatus (Squamata: Gekkonidae). - 
Zootaxa, 1493: 41-51. 


INGapy, M. Z. L. (2011): Habitat préférentiel d Uroplatus gigan- 
teus (Reptilia, Gekkonidae), dans le complexe de la Montagne 
d'Ambre, Nord de Madagascar. - Malagasy Nature, 5: 89-103. 


IUCN (2012): IUCN Red List Categories and Criteria: Version 3.1. 
- IUCN, Gland, Switzerland and Cambridge, UK. 


KATON, K. & D. M. STANDLEY (2013): MAFFT multiple sequence 
alignment software version 7: improvements in performance 
and usability. - Molecular Biology and Evolution, 30: 772- 
780. 


KUMAR, S., G. STECHER & K. TAMURA (2016): MEGA7: Molecular 
Evolutionary Genetics Analysis version 7.0 for bigger datasets. 
- Molecular Biology and Evolution, 33: 1870-1874. 


LIBRADO, P. & J. Rozas (2009): DnaSP vs: A software for compre- 
hensive analysis of DNA polymorphism data. - Bioinformat- 
ics, 25: 1451-1452. 

NUSSBAUM, R. A. & C. J. RAXwoRTHY (1994): A new species of 
Uroplatus Duméril (Reptilia: Squamata: Gekkonidae) from 
southern Madagascar. - Herpetologica, 50: 319-325. 


NUSSBAUM, R. A. & C. J. Raxwortuy (1995): New Uroplatus 
Dumeéril (Reptilia: Squamata: Gekkonidae) of the ebenaui- 
group from the Anosy Mountains of southern Madagascar. - 
Copeia, 1995(1): 118-124. 

Papu, J. M., A. MIRALLES, I. DE LA RIVA & M. Vences (2010): 
The integrative future of taxonomy. - Frontiers in Zoology, 7: 
article 16. 


PALUMBI, S., A. MARTIN, S. ROMANO, W. O. MCMILLAN, L. STICE 
& G. GRABOWSKI (1991): The simple fool’s guide to PCR. Ver- 
sion 2. — Honululu, Hawai. 


PUILLANDRE, N., S. BROUILLET & G. ACHAZ (2021): ASAP: as- 
semble species by automatic partitioning. - Molecular Ecolo- 
gy Resources, 21: 609-620. 


RATSOAVINA, F. M, E. E. Louis Jr., A. CROTTINI, R. D. RANDRI- 
ANIAINA, F GLAW & M. Vences (2011a): A new leaf tailed 
gecko species from northern Madagascar with a preliminary 
assessment of molecular and morphological variability in the 
Uroplatus ebenaui group. - Zootaxa, 3022: 39-57. 


RATSOAVINA, F, E Graw, N. RABIBISOA & N. A. RAKOTON- 
DRAZAFY (2011b): Uroplatus sikorae (amended version of 
2011 assessment). — The IUCN Red List of Threatened Species 
2011: e.T172820A6924215. https://dx.doi.org/10.2305/IUCN. 
UK.2011-2.RLTS.T172820A6924215.en, accessed 09.10.2022. 


RATSOAVINA, F. M., M. VENCES & E. E. Louis Jr. (2012): Phylo- 
geny and phylogeography of the Malagasy leaf-tailed geckos in 
the Uroplatus ebenaui group. - African Journal of Herpetolo- 
gy, 61 (2): 143-158. 

RatTsoavina, F. M., N. R. Raminosoa, E. E. Louis, JR., A. P. 
RASELIMANANA, F. GLAw & M. Vences (2013): An overview 
of Madagascars leaf tailed geckos (genus Uroplatus): species 
boundaries, candidate species and review of geographical dis- 
tribution based on molecular data. - Salamandra, 49: 115-148. 


RATSOAVINA, F. M., F. A. RANJANAHARISOA, F. GLAw, A. P. RASE- 
LIMANANA, A. MIRALLES & M. Vences (2015): A new leaf- 
tailed gecko of the Uroplatus ebenaui group (Squamata: Gek- 
konidae) from Madagascars central eastern rainforests. - 
Zootaxa, 4006: 143-160. 


RATSOAVINA, F. M., P. S. GEHRING, M. D. ScHERz, D. R. VIEITES, 
F. GLaw & M. VENCEs (2017): Two new species of leaf-tailed 
geckos (Uroplatus) from the Tsaratanana mountain massif in 
northern Madagascar. - Zootaxa, 4347: 446-464. 


RATSOAVINA, F. M., A. P. RASELIMANANA, M. D SCHERZ, A. RA- 
KOTOARISON, J. H. RAZAFINDRAIBE, F. GLAW & M. VENCES 
(2019a): Finaritra! A splendid new leaf-tailed gecko (Uropla- 
tus) species from Marojejy National Park in north-eastern 
Madagascar. - Zootaxa, 4545: 563-577. 


RATSOAVINA, F. M., M. D. SCHERz, K. A. TOLLEY, A. P. RASELIMA- 
NANA, F. GLAW & M. Vences (2019b): A new species of Uro- 
platus (Gekkonidae) from Ankarana National Park, Madagas- 
car, of remarkably high genetic divergence. - Zootaxa, 4683: 
84-96. 


177 


PHILIP-SEBASTIAN GEHRING et al. 


RATSOAVINA, F. M., F. Graw, A P. RASELIMANANA, A. RAKO- 
TOARISON, D. R. VIEITES, O. HAWLITSCHEK, M. VENCES & M. 
D. SCHERZ (2020a): Towards completion of the species inven- 
tory of small-sized leaf-tailed geckos: two new species of Uro- 
platus from northern Madagascar. - Zootaxa, 4895: 251-271. 


RATSOAVINA, F, F. Graw, N. RABIBISOA & N. A. RAKOTON- 
DRAZAFY (2020b): Uroplatus sameiti (amended version of 2011 
assessment). - The IUCN Red List of Threatened Species 2020: 
e. T193502A177286005. https://dx.doi.org/10.2305/IUCN.UK. 
20203.RLTS.T193502A177286005.en, accessed 09.10.2022. 


RAXWORTHY, C. J. & R. A. NUSSBAUM (2000): Extinction, and ex- 
tinction vulnerability of amphibians and reptiles in Madagas- 
car. - Amphibian and Reptile Conservation, 2: 15-23. 


RAXWORTHY, C. J., R. G. PEARSON, B. M. Zimxkus, S. REDDY, A. 
J. DEO, R. A. NUSSBAUM & C. M. INGRAM (2008): Continental 
speciation in the tropics: contrasting biogeographic patterns 
of divergence in the Uroplatus leaf-tailed gecko radiation of 
Madagascar. - Journal of Zoology, 275: 423-440. 


RONQUIST, F, M. TESLENKO, P. VAN DER MARK, D. L. AYRES, A. 
DARLING, S. HOHnAa, B. LARGET, L. Liu, M. A. SUCHARD & J. P. 
HUELSENBECK (2012): MrBayes 3.2: Efficient Bayesian phylo- 
genetic inference and model choice across a large model space. 
- Systematic Biology, 61: 539-542. 


SALZBURGER, W., G. B. EWING & A. VON HAESELER (2011): The 
performance of phylogenetic algorithms in estimating hap- 
lotype genealogies with migration. - Molecular Ecology, 20: 
1952-1963. 

SCHERZ, M. D., F. GLaw, C. R. HUTTER, M. C. BLETZ, A RAKO- 
TOARISON, J. KÖHLER & M. VENCES (2019): Species complexes 
and the importance of Data Deficient classification in Red List 
assessments: the case of Hylobatrachus frogs. - PLoS One, 14: 
€0219437. 

SCHNEIDER, J. G. (1792): Amphibiorum physiologiae specimen 
alterum, historiam et species generis Stellionum seu Gecko- 
num sistens. - Typis Chr. Ludov. Frid. Apitzii. Commissum 
Iohanni Andr. Kunze, Traieti ad Viadrum, Frankfurt a.d. 
Oder. 


SCHÖNECKER, P. (2008): Geckos Madagaskars, der Seychellen, 
Komoren und Maskarenen. - Terralog 12, Frankfurt/Main. 


SHEN, X. X., D. LIANG & P. ZHANG (2012): The development of 
three long universal nuclear protein-coding locus markers and 
their application to osteichthyan phylogenetics with nested 
PCR. - PLoS One, 7: e39256. 


STEPHENS, M., N. J. SMITH & P. DONNELLY (2001): A new sta- 
tistical method for haplotype reconstruction from population 
data. - American Journal of Human Genetics, 68: 978-989. 

SVATEK, S. & S. vAN DUIN (2002): Plattschwanzgeckos. Die Gat- 
tung Uroplatus. - Brähmer Verlag. 

Topp, M. (2011): Trade in Malagasy reptiles and amphibians in 
Thailand. - TRAFFIC Southeast Asia, Petaling Jaya, Selangor, 
Malaysia. 

UNEP-WCMC (2015): Review of selected species subject to long- 
standing import suspensions. Part I: Africa. - UNEP-WCMC, 
Cambridge. 

VENCES, M., A. Lima, A. MIRALLES & F. Graw (2014): DNA bar- 
coding assessment of genetic variation in two widespread 
skinks from Madagascar, Trachylepis elegans and T. graven- 
horstii (Squamata: Scincidae). - Zootaxa, 3755: 477-484. 

VENCES, M., A. MIRALLES, S. BROUILLET, J. Ducassr, A. FEDO- 
SOV, V. KHARCHEY, I. KosTapINov, S. KUMARI, S. PATMANIDIS, 


178 


M. D. SCHERZ, N. PUILLANDRE & S. S. RENNER (2021): iTaxo- 
Tools 0.1: Kickstarting a specimen-based software toolkit for 
taxonomists. — Megataxa, 6: 77-92. 


VENCES, M., S. PATMANIDIS, V. KHARCHEV & S. S. RENNER 
(2022a): Concatenator, a user-friendly program to concatenate 
DNA sequences, implementing graphical user interfaces for 
MAFFT and FastTree. - Bioinformatics Advances, 2. 


VENCES, M., M. MULTZSCH, S. GIPPNER, A. MIRALLES, A. CROT- 
TINI, P.-S. GEHRING, A. RAKOTOARISON, F. M. RATSOAVINA, 
F. GLAw & M. D. ScHERz (2022b): Integrative revision of the 
Lygodactylus madagascariensis group reveals an unexpected 
diversity of little brown geckos in Madagascars rainforest. - 
Zootaxa, 5179: 1-61. 


WILMÉ, L., S. M. GOODMAN & J. U. GANZHORN (2006): Biogeo- 
graphic evolution of Madagascar’s microendemic biota. - Sci- 
ence, 312: 1063-1065. 


WOLLENBERG, K. C., R. K. B. JENKINS, R. RANDRIANAVELONA, R. 
RAMPILAMANANA, M. RALISATA, A. RAMANANDRAIBE, O. R. 
RAVOAHANGIMALALA & M. VENCES (2011): On the shoulders 
of lemurs: pinpointing the ecotouristic potential of Madagas- 
cars unique herpetofauna. — Journal of Ecotourism, 10: 101- 
117. 


Supplementary data 


The following data are available online: 


Supplementary Table S1. Summary of locality information for 
samples of the Uroplatus sikorae complex confirmed by molecular 
data. 


